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The simplest non-trivial model of chaotic Bohmian dynamics is identified. We argue that its most 
important features can be observed in more complex models, above all, the presumable mechanism 
of the appearance of chaos in the Bohmian-type dynamical systems. 



I. INTRODUCTION 

t-H ■ 

In the de Broglie-Bohm approach to quantum phenomena J|, ||, ||| particles possess highly non-classical but well- 
. defined trajectories. They are derivable from the guidance equation: 
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where ^(r, t) = i?(r, t) exp [(i/h)S(r, t)] is a generator of the velocity field. 

Thus, so-called quantum chaos can be studied via Eq. (Q) since trajectories naturally exist in the de Broglie-Bohm 
mechanics which was proved to be completely equivalent to the standard Copenhagen version of quantum mechanics. 
The velocity field ((jjl preserves its definition also within the hydrodynamical formulation of quantum theory. 
' Solutions of Eq. ([j]) can be very complicated as it was first believed by Bohm himself Since the equations can 
be nonlinear ones, Diirr et al concluded that . . . there is nothing in Bohmian mechanics which would preclude 
sensitive dependence on initial conditions of . . . Bohmian orbits . . . and hence positive Lyapunov exponents. The idea 
of looking for chaos in Bohmian trajectories, and thus in quantum theory, was also mentioned in j|, |(J and some time 
later it was actually fulfilled j?J. Shortly after that a number of authors ||, [| |l(], |ll|, [l2|, [H| Q found chaotic orbits 
! | for various models and entailed wave functions. 

A generic feature of Bohmian mechanics is that the phase space volumes are not conserved by the flow, i.e.V • v 
does not generally vanish. It is expected however, that the volumes are conserved on the average, which means 
vanishing of limr_ >00 (l/T') f Q V-v dr. Very recently we have been successful [|l5| in constructing, within the Bohmian 
mechanics, a model of a Hamiltonian system (V • v = 0), perturbed by disturbance periodic in time, leading to chaotic 
£H . solutions for some parameters. 

Most studies on quantum chaology mentioned here g, [u], |ll|, [jj], |l3[ |l^| were inspired by Parmenter and 
Valentine's work Q|, where the system of two non-autonomous equations (or three autonomous), generated from 
the 2D anisotropic harmonic oscillator wave functions, was integrated. This work opened a number of interesting 
questions, among them: (i) what are the necessary conditions for the causal chaotic orbits to exist, especially, whether 
a model-system of equations playing the role of a limit case for other causal models can be identified, (ii) whether the 
known properties of volume-conserving systems with a time-dependent perturbation can be observed in the Bohmian 
dynamics as well. 

The problems we have just distinguished, though subjectively chosen, are however important for the causal dynamics 
itself and for the theory of dynamical systems in general. So far, they have not been fully answered. The reason is 
that the systems based on Eq. (Q) are much more difficult to deal with than any other conservative or dissipative 
system studied so far. In this paper we shall therefore try to identify a model being not only the simplest one in the 
causal dynamics but also revealing the way the chaotic Bohmian orbits are created. In what follows, the answer to 
the point (i) will be given and we shall also contribute to the point (ii) above. 



II. ESSENTIAL EQUATIONS 

A particularly interesting problem connected with the point (i) above is the following: what is the simplest form of 
the wave function V>(r, t) in Eq. (|l|) that still generates chaotic solutions. After having conducted numerical experiments 
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authors of the work |Q concluded that it is necessary for ip(r,t) to be a superposition of at least three stationary 
one-particle states and at least one pair of the states must have mutually incommensurate energy eigenstates. The 
wave function of the two-dimensional anisotropic harmonic oscillator used for the study generated a set of two non- 
autonomous equations with three control parameters. 

Using the same eigenstates we have recently shown ||l5| that a linear combination of only two stationary states is 
sufficient to obtain chaotic causal trajectories if one of the states is at least double degenerate. The dynamical system 
derived in such a way had two control parameters and when they were put equal to each other the system was proved 
to be completely integrable. 

We are going to show that with two properly chosen stationary states dynamical systems can be generated with one 
or two or four control parameters leading in each case to chaotic behaviour. To this end let us assume (h = 1, m = 
1, w = l) 

il>(x,y,t) = i(> n (x)ip n (y)e- iElt + [a^ k {x)^ n {y) + iaMx)^ k {y)] e"^ 2 ', (2) 

where tpj are solutions of the stationary ID Schrodinger equation and E k are eigenenergies of the 2D problems. The 
second stationary state is obviously double degenerate and ao and a\ are arbitrary real constants. From the definition 
of the velocity field given in Eq. ([[]) we have 

a (p'{x) [sin(ei) + anp(y)] 



(cos(ei) + a (p(xjf + (sin(et) + a^y)) 2 ' 
aif'iv) [cos(ei) + a tp(x)] 



(cos(et) + a ip(x)) 2 + (sin(et) + ai^(y)) 2 ' 



where 

V{Z) = P¥\, e = E 2 - El (4) 

and ip' means the first derivative with respect to z. 

From its construction the set (||) has two free parameters ao and oi. Nevertheless, their number can grow up to four 
or more, where E\ and E% are fixed, when some ideas of the so-called supersymmetric quantum mechanics are used. 
This method is based on using supersymmetry transformations to derive potentials V isospectral to a given one, say 
V, i.e. with exactly the same eigenvalue spectrum. For details of the procedure we refer the reader to literature (see 
e.g.[|l6|]). Here we are only going to present some formulae ready to use. 

Let {ip n (x)} and E n be respectively eigenstates and eigenenergies of a Hamiltonian with potential V{x). Define a 
function I[x) = J_ \ipo(z)\ 2 dz, where tpo{z) is the ground (nodeless) state. Then, for properly chosen values of the 

real parameter A, we can generate jl6| a new potential V{x; A) = V{x) — (d 2 /d.x 2 ) ln[I(x) -I- A] with (7 = 0,1,2... less 
bound states than V(x) has. Now the eigenfunctions ip{x; A) corresponding to V(x; A) in the strictly isospectral case 
of q — with E n = E n have the form 

7 , „ , , s 1 ( 1 A ( d ip'o(x)\ , 

i, n+1 (x; A) = i, n+1 (x) + - ^^-—^ j ^- - j (5) 

where n — 0, 1, 2 . . . and A>0orA<— 1 and the prim denotes the first order derivative. Now, the normalized wave 
function ip of the ground state reads 



teA) = ^±^o(*). (6) 
l[x) + A 

In the limit of A — > oo the "new" functions tp and potentials V reduce to the "old" ones. 

The procedure sketched above can be further generalized and one can construct an infinite number of isospectral 
families V(x; Ai, A2, ■ ■ ■ ) having identical bound-state energies. 

Since the velocity field in Eq. ([[]) is a functional of wave function we can modify properties of the dynamical system 
as in Eqs. (||), based on two stationary states, by increasing the number of control parameters. In the simplest case, 
instead of Eq. (||), we can now propose 

i>\n(x,y,t) =tp n (x; X)ip n (y; n)e~ lElt + 

r . i (7) 
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where /i>Oor/i<— 1. With the replacement ip — > we have from Eq. ([j]): 

. _ a, (p' x {x) [sin(ei) + a^^y)] 



(cos(ei) + a <p x (x) f + (sin(ei) + a^^y)) 2 ' 
ai^(y) [cos(ef) + a (p\(x)] 



(cos(et) + a <f\(x)) 2 + (sin(rf) + a^^y)) 2 ' 



where 

^ fc (x; A) „ ij>k(y,n) 
^aW = ? , , <P»{V) = -f-, r- (9) 

Of course as A — » oo and /x — * oo the set of Eqs. (||) reduces to that given in Eqs. 

Since for a single stationary state or for a linear combination of two non-degenerate stationary states dynamical 
systems resulting from Eq. (|l|) can be easily proved to be non-chaotic, we conclude that equations (^j or ffi) are the 
simplest non-trivial systems within the Bohmian mechanics with possibly chaotic dynamics. At this point it should be 
mentioned that chaotic behaviour can be generated from Eq. ([[]) even for just one stationary state ||l8|], it is necessary, 
however, to deal with two particles. 



III. THE SIMPLEST MODEL OF CHAOTIC CAUSAL DYNAMICS 



We are now ready to identify the model that could play the role of a limiting case of more complex models of 
Bohmian dynamics and still revealing their non-trivial properties. 

To this end let us observe that Eqs. (§) and (||) have the same form and the latter were derived to show a possibility 
of introducing additional control parameters to the system under consideration. Thus, for the time being, we may 
restrict ourselves to Eqs. (||). 

In the simplest non-trivial case we can choose 

<p{z) = Cz, e = 1 (10) 

in Eqs. (0), which follows from using two lowest states of ID harmonic oscillator, i.e., n = 0, k = 1 and ipo(z) — 
Aexp(— (l/2)z 2 ), ipi(z) — Bzexp(— (1/2) z 2 ). Due to the relation C = B/A, for normalized states C = \/2. Since 
the constant can be absorbed by ao and a\ we take for simplicity C = 1. Now ip(x,y,t) in Eqs. (||) is the state 
representing 2D isotropic oscillator with eigenenergies E\ = 1 and E<i = 2, 

In the case of ao = a\ = a, since the integral of the motion C(x,y,t) = M — a 2 \nM — 2a(xcost + t/sint), 
where M = (cost + ax) 2 + (sini + ay) 2 exists, Eqs. (||) do not generate chaotic solutions. Moreover, we also have 
dx/dx + dy/dy = which means conservation of the phase space "volume". It can be also observed that a function 
H(x,y,t) exists such that x = —dH/dy and y = dH/dx with H defined as H — (1/2) In M. 

When ao ^ a% Eqs. (||) are no more integrable and the very detailed discussion in [|l5) showed that the system of 
equations can be transformed into the form of a Hamiltonian autonomous system with a periodic non-Hamiltonian 
perturbation. Then, the method of Melnikov function can also be used to prove formally the existence of chaotic 
orbits. 

The model discussed here is, as yet, the only known model of Bohmian dynamics with the above properties. 
Presumably no similar case with such unique features can be found. To show that it can serve as a reference model 
for other models we have prepared a sequence of stroboscopic maps in Fig. [j] for a few pairs of ao ^ ai such that 
a = 1.0154 and oj = 0.985A with A = oo, 100, 10, 4, 2, 1. 

The first plot in Fig. [l] represents a circle which is then deformed into a curve formed of two loops crossing a single 
point. When ao = a\ both the homoclinic orbit of that shape and the hyperbolic point can be found formally as shown 
with full particulars in |L5|]. Here we show instead for A = 10 a thin stochastic layer appearing in the vicinity of the 
orbit spreading out in the area surrounding the point. The smaller the value of A is, the faster the homoclinic orbit 
breaks up and the fully developed chaotic trajectory can be attributed to the picture for A = 2. We have estimated 
the largest Lyapunow exponent for the trajectory with x(0) = 1.98, y(0) = and < t < 200000, At = 0.001, as 
Amaz = 0.06. 



IV. COMPARISON WITH MORE INTRICATE MODELS 



We shall now deviate from the above model to show that other more complicated models reveal similar behaviour. 
Two completely different ways of doing such modifications are briefly sketched below. 
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A. An isospectral modification 

We can destroy the rotational symmetry of the 2D harmonic oscillator potential of the model of Section 3 and 
simultaneously preserve the same eigenvalues spectrum of the new potential. To this end the functions ipo(z) and 
ipi(z), with z — x or y, are used in Eqs. (||) and (||) and will get ipo(x;X), ipi{x]X) and ipo(y; fi), ipi{y;i±). Finally, 



from Eqs. (||), we have: 



, , > _ I{x) + A ( e- x \ 

VxW ~ ^W+X) \ X + 2^{X + I(x)) J ' 

Mv) = -f=== [ y ■' 1 



where I(z) = 7r -1 / 2 J*^ e~" du. Once more the factor of a/2 is absorbed by the expansion coefficients in Eq. (|t]). 

The functions tpo(x; A) i^o(y; /J>) and ipa(x; A) 4>i{y; (J,) correspond to the same previous energies E\ — 1 and £»2 = 2, 
and hence again e = 1. The new potential has now the form 

V Xfi (x,y) =V x (x) + %{y) = 

I, o 2N 2 7 r- 1 /2 2 . e — 2 (A + /(.x))+^- 1 e- 2a;2 

- -(x +y) + 



2 V ~ ' " 7 ' (A + /(.x)) 2 (12) 

2^- 1 / 2 ;/e-^ (// + J(y)) + n^e-^ 2 
+ (M + %)) 2 

and in the limit of A — » oo and jj, — > cx> it reduces to the usual 2D oscillator potential. For the sake of comparison, 
its ID image is presented in Fig. ||[ The lower the value of A is, the more V\{x) deviates from its partner potential 
V(x) = (l/2)a: 2 . Similarly, <p\(x) in Eqs. ( |TT| ) deviates from <p{z) in Eq. (|Tc|). Thus, the model of Section 3 is a 
limiting case of the one introduced above for asymptotic values of A and /z. Only in this limit with ao = ai the latter 
is obviously integrable. In other cases, for any combination of the four parameters ao, ai, A, fj,, a chaotic solution can 
always be found. The components of the velocity field v(x,y, t) are now determined by Eqs. (||) and (|i~l|). Formal 
calculations show that V • v is now changing with time and tends to zero only in the above limit. Nevertheless, the 
scenario of approaching a chaotic orbit seems to be very similar to that observed for the model given by Eqs. ( [Tc| ) 
and (||). This is clearly visible in Fig. |[ 

The first two plots are about the same as in Fig. |l| and for smaller values of a the stochastic layers appear to be a 
little more pronounced than those in the corresponding pictures in Fig. [j]. We have expected this as a consequence of 
the lack of the rotational symmetry of the isospectral Hamiltonian (cf . Fig. |^) . 



B. Particle in a restricted space 

The second modification of our Eqs. ( fl"o|) and (||) uses the simplest model with a non-linear eigenvalue spectrum. 
It is a square well with x € [0, ir], y € [0,7r] and the normalized wave functions ip n = a/ (2/7t) sin (nz) and energies 
E n = (l/2)n 2 , n = 1, 2 . . . . Then, for n = 1 and k = 2 in Eq. @ we use in Eqs. (§) 

3 

ip(z) = 2cosz, e = —. (13) 

In our numerical calculations the factor of 2 is included in the constants ao and a\ in Eqs. (||). Now, Eqs. @ are 
integrated for a number of values of ao = a± = a and the results are presented in Fig. |[ The first member (|a| >> 1) 
of the family of pictures is not a circle but a closed curve perfectly approximated by the relation A = sinx • siny, 
where < A < 1. For smaller values of the parameter \a\ we can again observe the formation of a homoclinic orbit 
with the shape as in Fig. | for a = —10. When decreasing \a\ the orbit breaks up and the stochastic layer becomes 
well pronounced. 

The dynamics here is very similar to that in our model of Section 3 and this again suggests considering it as the 
simplest non-trivial model of chaotic Bohmian dynamics. 
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V. CONCLUSION 

In this article we have discussed the set of equations that could serve as the reference one for chaotic causal 
trajectories and revealing the way in which the causal orbits become chaotic ones. This work is the first attempt to 
contribute to this difficult problem. The results presented here seem to suggest that the way from a regular to chaotic 
behaviour leads via the formation and break-up of a homoclinic orbit. Unfortunately, the formal derivation of the 
orbit and of the critical points is possible only for the model of Section 3 above and for the details we refer the reader 
to |l5[ . In other cases we are able to do that only numerically. To emphasize this, we have prepared sequences of 
pictures in Figs. ^ and || clearly showing their close similarity to those in Fig. [[]. It is indeed common for the way the 
chaotic behaviour appears in Fig. |l| to be essentially independent of the model that is used to generate it. Results of 
our paper thus bridge the gap between the theory of dynamical systems and the causal quantum trajectories. 

We can consider the sequence of pictures in Fig. [j] as a possible way leading to the appearance of chaotic orbits 
in the Bohmian, or equivalently, in the hydrodynamical formulation of quantum mechanics. Quite strong support of 
the conjecture is in our opinion convincingly manifested in the series of plots in the above figures. What we have 
discussed here has not been observed so far since the models under consideration were too complex and the properties 
of trajectories depicted here were masked in their very complicated dynamics. 

The proof that some orbits in Fig. |[ and hence in Figs. || and [I] as well, are chaotic ones, follows directly from the 
results of our recent paper That is why we have calculated the largest Lyapunow exponent just for one orbit 

with A — 2 to show what the typical order of magnitude for the quantity and models under consideration is. 

Looking for the simplest yet non-trivial dynamics systems, both conservative and dissipative ones, has attracted 
much interest for a number of years. Representative studies on the subject together with the lists of such models can 
be found in fl~o| , [2(J. Our work is a contribution in this field for a special class of dynamical systems generated by 
the guidance equation ([]]) of the de Broglie-Bohm quantum mechanics. The model based on Eqs. (||) and ( |i~0| ) can be 
considered as the simplest non-trivial one for the reasons outlined throughout the paper. The name is additionally 
justified by the fact that the wave function we have used has only one moving node. Our numerical experiments with 
a number of other models lead to the conjecture that for the systems of two non-autonomous equations of the form 
of Eq. (jl]), existence of at least one moving node seems to be a necessary condition for chaotic orbits to exist. In our 
model the node moves along a circle in the integrable case and along an ellipse as the system is chaotic. 

We should emphasize at this point that a huge number of dynamical systems can be generated with help of Eq. 
(|l|) since they are a functional of the used wave function. The resulting systems of two non-autonomous equations, 
necessary by the well-known Poincare-Bendixson theorem for chaotic behaviour to appear, are not similar to any of 
systems studied so far. Our results show, however, that for the particular choice of the velocity field generator, i.e. the 
wave function, we can observe mechanism of obtaining chaotic orbits similar to the one known for the systems forced 
by the disturbance periodic in time. To find such a behaviour for "the quantum trajectories" of Eq. (|l]) we were 
constrained to a specific choice of functions without a particular physical significance. The similarity of Bohmian 
systems to the well-known dynamical ones is presumably still preserved for physically important wave functions. This 
statement is, for now, an open question. 
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FIG. 1: Stroboscopic maps of the period 2ir for the model described by Eqs. (|3|) and (fL0|). In each case the time interval 
< t < 50000 with the time step of At = 0.001. 
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FIG. 3: Stroboscopic maps of the period 2n for the isospectral model described by Eqs. (g) and ( |ll| ) with do = a\ — a and 
A = fi = 20. The time interval and the time step are respectively < t < 50000 and At — 0.005. 
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FIG. 4: Stroboscopic maps of the period (4/3)7r for an infinite square well model described by Eqs. (|3j) and (fL3|). The time 
interval is in each case < t < 50000 and At = 0.001. 



